Quantum chaos and operator fidelity metric 
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We show that the recently introduced operator fidehty metric provides a natural tool to investigate 
the cross-over to quantum chaotic behaviour. This metric is an information-theoretic measure of 
the global stability of a unitary evolution against perturbations. We use random matrix theory 
arguments to conjecture that the operator fidelity metric can be used as an "order parameter" to 
discriminates phases with regular behaviour from quantum chaotic ones. A numerical study of the 
onset of chaotic in the Dicke model is given in order to support the conjecture 



Introduction- The information-theoretic task of distin- 
guishing quantum states can be formulated in terms of 
the differential-geometric notion of metric over the quan- 
tum state space both in the pure and mixed state case 
[l[ . This remarkable fact lies at the root of the recently 
proposed fidelity approach to quantum phase transitions 
(QPTs) 3- The intuition behind this approach is quite 
simple: at the transition points a slight change of the 
system driving parameters e.g., external fields, gives rise 
to a major structural change of the associated (thermal 
or ground) quantum state and this, in turn, should result 
in an enhanced statistical distinguishability. In view of 
the above mentioned connection one expects detectable 
signature in the metric function e.g., singularities in the 
thermodynamical limit. This intuition has been con- 
firmed by the study of a host of different many-body 
systems, at zero and finite temperature [1, 

In Refs |5j the ground state approach has been ex- 
tended to the operator level. In order to do that one has 
simply to notice that finite-time quantum evolutions over 
the Hilbert space Ti. correspond to unitary operators that 
belong to a space i.e., the space of linear operators C{TL), 
that can be turned itself into a Hilbert one in many dif- 
ferent ways. This allows one to introduce the notion of 
operator fidelity and the associated metric Q . 

The main goal of this paper is to attempt to use oper- 
ator fidelity to investigate another fascinating and some- 
what elusive phenomenon : quantum chaos [7|. In many 
years of theoretical efforts a variety of approaches, like 
random matrix theory quantum motion reversibility 
and stability up to very recent measures of phase- 
space growth rate have been proposed fll'l . All of these 
related approaches appear to be partially successful, yet 
a unified fully and satisfactory characterization of quan- 
tum complex behaviour is still missing. 

In the following we first review the relevant background 
on operator fidelity metric and we then carry out our 
analysis along two well separated though mutually sup- 
porting paths. In the first analytical part we explore 
the consequences of the physical intuition that is behind 
the approach that we advocate: the quantum chaotic 
evolutions can be characterized as those who have a re- 
silience against random perturbations. This resilience is 
here quantified by the information-theoretic distance be- 



tween perturbed and unperturbed evolutions. We will 
make use of random matrix theory [8,] arguments to for- 
mulate and support our main conjecture on operator met- 
ric (in)sensitivity of (chaotic) regular system to random 
perturbations. The second part of our paper offers an 
implementation these ideas in a particular system. More 
precisely, we numerically analyze the cross-over to quan- 
tum chaos in the Dicke model [H, E^l ^tnd we show how 
the operator metric approach can be used as a tool to 
identify and characterize this cross-over. We will con- 
clude by a summary and an outlook. 

Information- Theoretic metric over manifolds of 
unitaries- Let us start by providing the setup and the 
basic facts about operator fidelity metric found in 0]. 
Let 7i be a quantum state-space and p a density matrix 
over it (p > 0, trp = 1). One can define the following 
hermitean product over CiTi.): {X,Y)p :— tr(pX^Y). If 
p > then this scalar product is non-degenerate and 
\\X\\p := ^y (X, X)p defines a norm over C{H). The 
operator fidelity X,Y Cz C{Ti.) is then given by 



Tp{X,Y) ■.= \{X,Y)p\ 



(1) 



This quantity has well-defined operational meaning in 
terms of statistical distinguishability (the higher T the 
lower this latter) in terms of bi-partite quantum states 
, interferometric schemes associated with X, Y and p 
and decoherence rates [!'] . The case we will be concerned 
with is when X ^ Ux ■■= e-'*"\ Y = Uy{X ^ A') we 
will also assume that [p, H\] = 0. 

In passing we would like to note that ^ has relevance 
to non-equilibrium dynamics e.g., quenches where one 
goes suddenly from H\ to H'^^. Indeed H]) is the mod- 
ulus of characteristic function of the work distribution 
function [13] TiUxit),Ux>it)) = | / dwP(w)e*'^*| where 
Pi^) ■=Em,nPn,n\in\m')\^S{co-E^ + E'^,) (£;„ and |n) 
and the corresponding primed quantities denote eigenval- 
ues and eigenvectors of Hx')- Moreover for p„^„ — Sn,o 



one obtains a Loschmidt echo |(0|e' 



itH, 



pointing out 



to relevance of H]) to dephasing experiments ^, • 

In this paper, following the differential-geometric spirit 
of [2] we are going to consider the operator fidelity ^ be- 
tween infinitesimally different unitaries i.e.. A' = A -I- dX. 
The proof of the following proposition is just a direct cal- 
culation analog to the one performed at the state-space 
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level Let {U\} C U{T-L) be a smooth family of uni- 
taries over Ti parametrized by elements A of a manifold 
M. One finds Tp{Ux,Uy,+sx) = 1 - l2xp{\){U\U') 
where U' — dU / dX and 

Xp{X)iX,Y) := {X,Y)p-{X,U^)p(U^,Y)p (2) 

We see that the leading term in the expansion of ([T|) 
defines a quadratic form over the tangent space of the 
projective >C(7i) at U\. For full rank p that quadratic 
form is a metric. This operator metric (also referred to 
as operator fidelity susceptibility in 0]) is the central 
tool of this paper. The intuitive meaning of x is quite 
simple: the larger the operator metric at U\ the greater 
the degree of statistical distinguishability between the 
quantum evolutions associated with two slightly different 
set of control parameters A and \ + d\. This fact has been 
used in (Q) to use ^ to study quantum criticality. 

For the developments of this paper it is convenient 
to make use of the superoperator formalism common in 
Quantum Information Science. We define over C(Ti.) the 
superoperator Lh ■— [H^']. This is nothing but the 
generator of the Heisenberg evolution and it is easy to 
see that, if |n)'s and the E^s denote the eigenvectors 
and eigenvalues of H{X), En — E^ and \n){m\ are the 
eigenvalues and eigenoperators of Lh respectively. It 

^ 1/2 

follows that \'^n,m) = Pm,m\n){m\ it is an orthonormal 
basis of i.e., {i'n,ni\'^p,q) = 5n,p5m..,q- If wc define 

P = ^„ |^n,n)(^n.ri| as the projcctiou over the kernel 
of Lff, Q = i - P, and 5t[x) := [fiu\{xt / 2) / {x / 2)f , one 
finds that the operator metric is given by the sum of two 
distinct terns: xp = Xp^^ + X'p'^ i where 

X« = mLH)Q\H')\\l; ^ = \\P\H')\\l-\{t,H')p\^ 

(3) 

Operator fidelity and random matrix theory- From the 
last Eqs. we see that, in general Xp depends on t, H\ 
and on its derivative H' = dH\/dX. We would like to 
start our analysis of the chaos-related properties of Xp 
by introducing a ^-dependent quantity Xp that, contains 
information just on the Hamiltonian H\ . A possible way 
to achieve this goal is i) replace H' in the operator met- 
ric by a perturbation V that is assumed to be drawn 
by a Gaussian ensemble, with measure 2?[V^], of random 
matrices^, ii) take Xp{t,H) :— J 'D[V]xp{t, H,V) as the 
average Xp over the ensemble of Vs. 

More precisely, we will assume that the the pertur- 
bation matrix elements Vn.m are gaussian random vari- 
ables satisfying the relations / 'D[V]Vn.mVp^q 5n,p5m,q 
The key step is the following fact: If Rp denotes the 
superoperator X — > Xp one has / I?[y]|V")(y| ~ Rp 

[l^ . Now, by using this result, and writing Xp^"* = 

'Iv[\V){V\5t{LH)Q] and £^ = Tr [|^^)(y|(P - |l)(a|)] , 
one finds 

X^p{t,H) ^ Tv[Rp5t{LH)Q] , xf\t,H) ^ t^l - trp^) 



From these relations we clearly see that the two terms of 
the operator fidelity metric behave quite differently upon 
averaging over the perturbation V : on the one hand Xp 
looses any direct connection with the Hamiltonian iJ, 
just the purity of the state p is returned, on the other 
hand the averaged Xp^'' still shows an highly non-trivial 
dependence on the spectral features oi H . In the rest of 
the paper we are going to focus on this latter term. 

Let us make the content of the x'p^ explicit by re- 
sorting once again to the eigen-operator basis | of 

Lh, from Q we find X}t\t,H) ^ Y^n^m Pn,n5t{^n,ra) 

A„^„i :— En — Em- This equation shows that the pertur- 
bation averaged x'"'^' depends on the the of level spacing 
distribution of H. In particular - in view of the prop- 
erty YiTat^oot^^5t{x) — 2t:6(x) - for sufficiently large t 
we expect contributions from small A„^m i-e., by almost 
crossing levels, to dominate the operator metric. The 
key observation at this point is that one of the possible 
characterization of the presence of chaos in quantum sys- 
tem stems from the level spacing analysis ITI. Il3l|. In this 
context, the properly normalized spectrum |15j is ana- 
lyzed in terms of the spacing between consecutive energy 
levels 5'„ — En+i — En- If P{S) denotes the probability 
that two nearest-neighbouring levels have an energy dif- 
ference S then the behaviour P{S) for S* — + encodes 
the main features distinguishing a chaotic system from 
a regular one. Indeed, the distributions characterizing 
the regular spectra are Poissonian, allowing for a non 
vanishing probability to have S = Q spacing i.e., level 
crossings. In the chaotic case one has Wigner-like dis- 
tributions Pw{S) oc 5"'exp— S*^; their small spacing be- 
haviour, in view of symmetry, is dominated by S'^ , with 
V > Q. This entails for the phenomenon of level repulsion 
i.e., suppressed level crossings. 

The above remarks together put us now in the position 
to formulate the main result of this paper in the form of 
a Conjecture: for sufficiently large t the operator met- 
ric Xp^^ is finite for regular quantum Hamiltonian and 
(almost) vanishing for quantum chaotic ones. 

To provide further support to this conjecture let us 
consider the average behavior of Xp^^ with respect H i.e., 
X'p\t) := / V[H]x'p\t,H). Since x'p\t, H) is a spectral 
function integration over the eigenvectors of H can be 
readily performed and one is left with ^] 

xi^\t)^ f e-^"^" H dEndEm{En-Emrxi^\t,H) 

n<m 

The exponent ly tells apart ensembles of regular, 1^ = 0, 
from chaotic, v > Hamiltonians. From the last equa- 
tion, and by resorting to the explicit representation Xp^'' 
above we conclude liuit^oo Xp^\t)/t ^ for regular H 
and Ihnt ^ oo x^p\t) /t « for chaotic H. Assuming typ- 
icality i.e., typical and average behavior coincide, in H 
these Eqs are nothing but a formulation of the Conjecture 
above. 
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Of course none of the above arguments is rigorous, for 
as they depend on a set of unproven assumptions e.g., 
RMT ensembles physical relevance, typicality, moreover 
they somehow neglect the potential sensitivity on p and 
finally rely on a "naive" large t limit. Nevertheless the 
conceptual content of the Conjecture is quite compelling 
and intuitive at the same time: typical regular Hamil- 
tonians have a much higher susceptibility against ran- 
dom perturbations than chaotic ones. In other words, a 
random perturbation to a regular quantum evolution re- 
sults in new evolutions that has typically a much greater 
information-theoretic distance from the original one than 
if the same perturbation were applied to a quantum 
chaotic evolution. 

Operator fidelity and transition to chaos: the Dicke 
model- In this section we complement the above theoret- 
ical analysis of the operator fidelity with a numerically 
study. We will use the fidelity metric as a tool to study 
the transition from regular to a chaotic regime in the 
Dicke model [l^, ■ This model consists of a set of N 
identical spin 1/2 atoms with atomic level splitting ujq 
placed in an ideal cavity that collectively interact with a 
single bosonic field described by the operators {o, a'''} = 1 
and characterized by the frequency uj = ojq (resonance). 
The Dicke Hamiltonian reads 

= WoJz + wa^a + ^= (a^ + a) (J+ + J_). (4) 

where = Yh=i '^i and J± = ± Jy = J2iLi '^f ^-re 
coUctive spin operators. The pseudo-spin length is fixed 
(j — N/2) and we have that the interaction is weighted 
by = 1/^/N. The Dicke model is integrable at 

finite sizes only if, for small A, one implements the rotat- 
ing wave approximation (RWA) i.e., if one neglects the 
terms J-i-, oJ-. In the thermodynamical limit (TDL) 
iV — > oo [l3 . [Tsj it is integrable and it exhibits a quan- 
tum phase transition in correspondence of A = 0.5. This 
critical value separates the "normal" phase (A < 0.5) and 
the "super-radiant" phase and the related phase transi- 
tion can be appropriately described in the context of the 
state fidelity approach 

Here we are interested to the finite size instance of the 
model that, in absence of the RWA is not integrable and 
exhibits a regular to chaos cross-over. Indeed, the only 
(known) symmetry property is described by the parity 
operator P = exp[i7rN], where N = a'^a+Jz+j is the op- 
erator that counts the number of total quanta present in 
the system. In order to characterize the quantum chaotic 
behaviour of Q one can then resort to the study of the 
properties of the spectrum relative to the odd (even) sub- 
spaces. Indeed, in T3] the authors showed that in cor- 
respondence of the value of A — 0.5, which is critical in 
the TDL, one can observe a transition from a regular re- 
gion (A < 0.5) characterized by Poisson-like level spacing 
distributions, to a chaotic region (A > 0.5) where the 
distributions are Wigner-like with v ^ 1. 

The actual system we used in our simulations is a set 
of = 20 atoms coupled with a bosonic bath. To make 
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FIG. 1: Plots of Xp^'(r,A) corresponding to T = 1 (a) and 
T = 10 (b). The state is p = Inset: plot of the relative 

entropy S{Px\\Pw) 



computations feasible we have to introduce a cut-off in 
the -infinite dimensional ~ bosonic state space. This cut- 
off has to be chosen in such a way that the bosonic system 
still operates as a bath for the atoms. We thus choose to 
include the first Ad — 128 bosonic energy levels; the total 
Hilbert space have thus dimension d — 4032. 

Let us first consider the level spacing distribution 
P\{S) of H{X). A characterization of P\{S) can be 
given by its statistical distinguishability with respect to 
the Wigner distribution P\y{S) {v — 1); this statisti- 
cal distance can be measured by the relative entropy 
S{Px\\Pw) ■■= E„^A(5„)log[FA(5„)/Fw(^„)] 0. In 
the inset of Fig. ((TJa) we see that in the chaotic chaotic 
region the relative entropy is uniform and small i.e., 
Px{S) ~ Pw{S), whereas in the regular one we have 
in general higher values and a more complex behavior 
since the degree of distinguishability varies with A. The 
transition point being roughly at A = 0.5. 

We now analyze the numerical results for x'/P (^i ^) = 
'En^mPn.n\{n\H'\n)\-'StiA^^^) whcrc H' = dxH =. 

\/2j (flt -|- a) (J+ -|- J_) implements the perturbation to 
H. We have first chosen p = 1/d i.e., scalar product in 
([1]) is Hilbert-Schmidt, this corresponds to an infinite- 
temperature thermal state over the truncated working 
space. In Fig. ([T]) and ^ we show Xp'' (t ~ T, A) (nor- 
malized to its maximal value max^ Xp^^(T, A)) for differ- 
ent T's. These plots clearly show that for growing values 
of T the behavior of Xp^^C^i changes significantly. For 
T > 1 the response of the system to the infinitesimal 
change of parameter A, as described by x'p\tj clearly 
exhibits two markedly different behaviour for A < 0.5 and 
A > 0.5. Indeed this behaviour is coherent with the con- 
jecture described in the previous paragraphs: the long 
time response of the system is characterized by the sen- 
sitivity of X^pHt, A) to the S* — *■ part of the spectrum. 
In particular, in the regular region, consistently with the 
relative entropy analysis, the behaviour is non uniform 
with A and the (non uniform) presence of level crossings 

give rise to a pronounced response of Xp^^ characterized 
by a sequence of high peaks. In contrast, in the chaotic 
region the level repulsion phenomenon is reflected in a 
response which is uniform and low for A > 0.5 
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FIG. 2: Plots of xi^HT,X) corresponding to T = 100 (a) and 
T = 1000 (b). The state is p = 1/d. 
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FIG. 3: Left (Right) plot of Xp^' (^i ^) corresponding to t — 
100{t = 1000). The state is p = exp (-/3H)/Tr[exp (-/3i?)] 
where P = 0.014. 



We finally complete our analysis by showing in Fig. 
([3]) how Xp^"* (i, A) behaves by choosing in ([1]) the ther- 
mal state p = exp (— /3i?)/Tr[exp (— with finite /?.. 
In principle the introduction of a finite temperature, by 
inducing a relative weight among the energy levels could 
spoil their chaotic vs regular distributions. Here we have 
chosen the inverse temperature /3 = 0.014; this corre- 
sponds to fixing the ratio between the thermal probabil- 



ities relative to the ground and the highest energy state 
ex-p[—(3{E,nin — Emax)] ~ 0.05. Morcovcr we checked 
that for this /? the spectral features of P{S) relevant 
to chaos are not washed out. Fig. ^ shows that the 
T = 100 and T = 1000 behaviour of x'p'' {t, A) is again 
markedly different in the two regions; the response of the 
system is enhanced by the presence of level crossings in 
the regular region, while it is suppressed in the chaotic re- 
gion. These results are again consistent with our conjec- 
ture and suggest that the operator fidelity metric might 
well work as an order parameter for quantum chaos in 
appropriately chosen temperature intervals. 

Conclusions- In this paper we have tackled the prob- 
lem of characterizing the chaotic properties of quantum 
systems by means of an information-geometrical tool: the 
operator fidelity metric. The results of our analysis are 
twofold. At a purely analytical level, by means of tech- 
niques drawn by random matrix theory, we have formu- 
lated and substantiated the conjecture that the operator 
fidelity metric may provide an order parameter for the 
regular to chaotic cross-over. The intuition being that 
a generic chaotic system, at variance with a regular sys- 
tem, is resilient with respect to random perturbations. 
We have then shown how the operator fidelity approach 
can be used as a tool to numerically identify the regular 
to chaotic cross-over in a relevant many-body system i.e., 
the Dicke model. We believe that the results obtained in 
this paper while preliminary, are promising and give rise 
to several interesting questions. Assessing the generality 
and efficiency of the methods we introduced and unveiling 
their relations with other approaches to quantum chaos 
e-g-, is the subject of ongoing investigation. 
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